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The phenomenon of current sheet canting in pulsed electromagnetic accelerators is the departure of the 
plasma sheet (that carries the current) from a plane that is perpendicular to the electrodes to one that 
is skewed, or tipped. Review of pulsed electromagnetic accelerator literature reveals that current sheet 
canting is a ubiquitous phenomenon - occurring in all of the standard accelerator geometries. Develop- 
ing an understanding of current sheet canting is important because it can detract from the propellant 
sweeping capabilities of current sheets and, hence, negatively impact the overall efficiency of pulsed elec- 
tromagnetic accelerators. In the present study, it is postulated that depletion of plasma near the anode, 
which results from axial density gradient induced diamagnetic drift, occurs during the early stages of the 
discharge, creating a density gradient normal to the anode, with a characteristic length on the order of the 
ion skin depth. Rapid penetration of the magnetic field through this region ensues, due to the Hall effect, 
leading to a canted current front ahead of the initial current conduction channel. In this model, once the 
current sheet reaches appreciable speeds, entrainment of stationary propellant replenishes plasma in the 
anode region, inhibiting further Hall-convective transport of the magnetic field; however, the previously 
established tilted current sheet remains at a fairly constant canting angle for the remainder of the dis- 
charge cycle, exerting a transverse J x B force which drives plasma toward the cathode and accumulates 
it there. This proposed sequence of events has been incorporated into a phenomenological model. The 
model predicts that canting can be reduced by using low atomic mass propellants with high propellant 
loading number density; the model results are shown to give qualitative agreement with experimentally 
measured canting angle mass dependence trends. 
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Nomenclature 

B - magnetic inductance [Gauss] 

J - current density [A/cm 2 ] 

u d - diamagnetic drift velocity [cm/s] 

v ba - branch/anode interface velocity [cm/s] 

v bt - branch/trunk interface velocity [cm/s] 

c - speed of light in vacuum [cm/s] 

e - electronic charge [statcoulomb] 

I sp - specific impulse [s] 

L - trunk to anode separation [cm] 

L c - characteristic density gradient length [cm] 
n e - electron number density [#/cm -3 ] 
nj - ion number density [#/cm~ 3 ] 
mi - ion atomic mass [g] 

R c - convective skin effect Reynolds number 

s - branch/trunk to branch/anode axial separation [cm] 

ha - Alfvdn speed [cm/s] 

v c - Hall-convective speed [cm/s] 

i>d - resistive diffusion speed [cm/s] 

a - system mass per unit power [kg/kW] 

Si - ion skin depth [cm] 

T) - plasma resistivity [11m] 

(9 - canting angle [degrees] 
cjp - plasma frequency [rad/s] 

1 Introduction 

P ulsed electromagnetic accelerators are devices which 
use intense bursts of electrical current (d(10 4 — 
10 fi ) A) to create high speed (0(1O 3 - 10 s ) m/s) jets 
of plasma. They find application as plasma sources in 
many basic plasma science experiments[l] as well as in 
a specific genre of electric space propulsion device called 
the pulsed plasma thruster (PPT)[2]. The present work 
is motivated by the desire to improve the performance 
of pulsed electromagnetic accelerators in the context of 
plasma propulsion. 

PPTs have the potential for fulfilling the attitude con- 
trol requirements on a satellite at greatly reduced mass 
and cost. They are also being considered for constellation 
maintenance for missions such as interferometric imaging 
of the Earth from space or deep space from an Earth orbit 
(c.f., Polzin et a/[3]). The benefits of PPTs are their sim- 
plicity, very small impulse bits for precise control of satel- 
lite motion, reliability, and high specific impulse. Two 
classifications of PPTs exist, corresponding to the form 
of propellant used: gas-fed (GFPPT) or ablative propel- 
lant (APPT). The gas-fed variety has the advantages of a 
“clean” exhaust plume and high specific impulse. The ab- 
lative version of the PPT uses a solid propellant, such as 
Teflon, to provide other advantages such as compactness 


and overall ease of system integration; however, plume 
contamination and lower specific impulse may limit the 
application of APPTs for some missions. The model pre- 
sented in the present study is relevant to the GFPPT, where 
electromagnetic acceleration is the dominant mode of ac- 
celeration; the results of this study may not be directly 
relevant to the APPT, where electrothermal acceleration 
may be significant[4]. 

2 Review of the Problem 

2.1 Definition of the Problem 

The phenomenon of current sheet canting is the departure 
of the current sheet from perpendicular attachment to the 
electrodes to a skewed, or tipped, attachment. It is best il- 
lustrated by an example. Figure 1 shows the evolution of 
a hypothetical and a real current sheet near the breech of a 
rectangular-geometry pulsed electromagnetic accelerator; 
outlines of the electrodes (the cathode is the bottom elec- 
trode) have been added for clarity. Ideally, we would like 
the current sheet to initiate at the breech, perpendicular to 
the electrodes, and remain so as it propagates axially (as 
illustrated on the left-hand side of the figure). In contrast, 
the right-hand side of the figure shows the experimentally 
observed evolution of a discharge^]. As expected, the 
current sheet is seen to initiate at the breech; however, as 
time progresses, the current sheet is seen to severely tilt, 
or cant, as it propagates. 

Canting creates off-axis components of thrust, which 
constitute a profile loss. Consider the t = 8 (is photograph 
on the right-hand side of Fig. 1 . Assuming that the mag- 
netic field is uniform behind the current sheet, the force 
on the top electrode (which is found by integrating the 
magnetic pressure on the electrode surface) will be greater 
than on the bottom electrode, because the top electrode 
has more surface area exposed to the magnetic field. This 
transverse force imbalance may result in an undesirable 
torque on a spacecraft which uses a PPT. Also, consid- 
ering the work done by the current sheet, it is clear that 
a canted current sheet will apply a force to the propellant 
transverse to the thrust axis and thus expend energy which 
is not converted into useful thrust. 

In addition to causing an off-axis component of thrust, 
current sheet canting may undermine the effective sweep- 
ing up of propellant as the current sheet propagates. The 
effect of canting may be to force the plasma entrained by 
the current sheet into the cathode where it stagnates and is 
then left behind. Indeed, canted current sheets may act, 
undesirably, like real snowplows - never accumulating 
but, rather, throwing their load to the side as they pass 
by. 
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t = t4 t = 1 0 [|is] 

Ideal Experiment 


Figure 1 : Comparison of ideal and experimentally observed current sheet propagation (from Markusic[5]). In the ideal 
case, it is implied that time tl <t2<t3<t4. The experimental photographs are from an argon (100 mTorr, uniform fill) 
discharge with peak current about 60 kA. Outlines of the electrodes have been added for clarity; the vertical rectangular 
element on the right-hand side of the pictures is a structural element that obstructed optical access. 


While direct studies of performance degradation due 
to current sheet canting are needed, the potentially ad- 
verse effects envisaged above provide reasonable justifi- 
cation for pursuing a theoretical study of the phenomenon. 
By developing an understanding of the physical processes 
which drive current sheet canting, we can develop pre- 
scriptions for how to reduce the effect, and ultimately 
provide guidance for the design of better pulsed plasma 
thrusters. 

2.2 Review of Past Research 

The literature from the early GFPPT researchers (see, for 
example, [7]-[15]) indicates that current sheet canting was 
a ubiquitous phenomenon - occurring in a variety of dif- 
ferent electrode geometries and experimental conditions. 
However, detailed treatment of the subject is limited, with 
most references being anecdotal in nature. A review of 


the earlier studies can be found in Markusic[5]; the con- 
clusions from that review can be summarized as follows: 

>• Current sheet canting always occurs in an orienta- 
tion such that the anode current attachment leads the 
cathode current attachment. 

>• Current sheets are always observed to cant with ni- 
trogen and argon propellant - irrespective of geom- 
etry (rectangular, coaxial, z-pinch, inverse z-pinch, 
and parallel rod), gas pressure, or current level. 

>■ Uncanted hydrogen and deuterium current sheets 
have been observed but only with specially prepared 
electrodes. 

>■ Ion current conduction is believed to play an impor- 
tant role in establishing the overall current pattern. 

>■ Current sheets appear to exhibit invariant features 
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(e. g. anode foot, cathode hook, and canting) over 
a wide range of geometries, propellent species, gas 
pressures, and current levels. 

>■ Conflicting experimental results exist with regard 
to the influence of polarity in coaxial accelerators, 
but variations in the construction (e.g., slotted elec- 
trodes) may account for the different observations. 

3 Model Description 

Below we present a phenomenological model for cant- 
ing that is motivated by the observations listed above as 
well as recent experimental measurements of current sheet 
canting angles for a variety of propellants[5]. 

3.1 Current Conduction Phases 

The life of a current sheet can be divided into three phases: 
starvation, branching, and snowplow phases. Figure 2 
schematically illustrates this evolution. Corresponding 
photographs from the experimental study[5] are shown at 
the top of the figure. The grey objects in the illustrations 
are meant to illustrate the spatial extents of the current 
sheet plasma. The solid lines behind the current sheet are 
meant to represent magnetic flux tubes (i.e., the amount 
of magnetic flux contained between consecutive pairs of 
lines is constant from frame-to-frame). 

The first illustration in Fig. 2 shows the initiation of 
the current sheet at the breech of the accelerator. In il- 
lustration 2 the current sheet is shown shortly after initi- 
ation. The plasma in the region near the anode becomes 
severely depleted due to mass motion of the fluid toward 
the cathode (for reasons to be explained later). This first 
phase is termed the “starvation phase”, borrowing termi- 
nology from a phenomenon observed in MPD thrusters 
called “anode starvation”[16] (which probably originates 
from an altogether different physical mechanism). 

Between illustrations 2 and 3 the plasma near the an- 
ode becomes so tenuous that it can no longer contain the 
magnetic field behind it. In a sense the current sheet 
can be thought of as a thin membrane (like a balloon) 
that contains a high pressure “magnetic fluid”. When the 
plasma near the the anode becomes sufficiently “thin”, the 
membrane quickly expands, or ruptures, allowing the con- 
tained magnetic flux to rapidly stream through. The ini- 
tial current channel (henceforth referred to as the “trunk”) 
becomes bifurcated along the anode, as the streaming 
magnetic flux (and associated surface current) propagate 
ahead of the original current attachment point; this inter- 
face forms a new conduction path for the current sheet 
which connects the trunk to the anode (this new conduc- 
tion path will be referred to as the “branch”). The newly 


formed branch and trunk form what was referred to by 
earlier researchers as the “anode foot”. Eventually the 
trunk becomes magnetically insulated (as the magnetic 
field wraps around the top of the trunk, transport from 
the top of the trunk to the anode is impeded by the trans- 
verse magnetic field) from the anode whereupon all of the 
current flows through the branch. 

The branch propagates, borrowing from shock-tube 
parlance, as a contact discontinuity, with the anode as one 
wall and the trunk as the other. The magnetic pressure 
drives the contact point down along the branch-trunk in- 
terface and forward along the branch-anode interface, re- 
sulting in a canted current sheet. Anode starvation occurs 
to a lesser extent in the branch because it continually prop- 
agates into a fresh supply of propellant (the branch-anode 
interface is replenished with propellant from the dynamic 
pressure associated with its substantial axial speed). On 
the other hand, the initial current sheet (the trunk), is un- 
able to avoid anode starvation because of its slow initial 
speed. As a final note on the branching phase, the mag- 
netic pressure between the branch and the trunk causes 
the trunk to deform into the hook-like structure reported 
in many studies; the hook is simply a vestige of the initial 
current sheet. 

The branching phase ends when the bottom of the 
branch reaches the cathode and the final phase (as shown 
in illustration 3), the snowplow phase, begins. The mag- 
netic pressure is uniformly distributed on the back face 
of the canted current sheet and the current sheet is suffi- 
ciently dense in all areas to prevent further field leakage. 
The current sheet thus remains at a fairly constant cant- 
ing angle during the remainder of its propagation. The tilt 
of the current sheet causes it to exert a cathode-directed 
component of the JxB force density on all of the propel- 
lant which is subsequently swept up. As a result, the pro- 
pellant is preferentially directed toward cathode, where it 
accumulates. This “mass-funnelling” to the cathode may 
cause elevated plasma pressure along the cathode - lead- 
ing to expansion of the propellant into the region behind 
the current sheet. A structure in the form of a plasma 
“bubble” behind the current sheet could form and, over 
time, grow large enough to span the entire gap between 
the anode and cathode and cause a “restrike”, which ef- 
fectively short circuits the initial current sheet. 

3.2 Terminal Canting Angle 

Current sheet canting results from a disparity in the axial 
propagation speed between the branch-anode and branch- 
trunk interfaces. The branch-anode interface moves faster 
and the current channel connecting the anode to the trunk 
forms a canted current sheet. 
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Figure 2: Illustration of a model showing current conduction phases which lead to a canted current sheet. 


The canting angle, 6(t), continuously evolves from 
zero to its terminal (final) value, Of, during the branch- 
ing phase. The severity of the canting depends on the 
relative magnitude the branch-anode and branch-trunk in- 
terface velocities. The terminal canting angle can be pre- 
dicted by calculating the trajectories of two points on the 
branch: the branch-trunk interface point (point BT) and 
the branch-anode interface point (point BA). The elements 
of the model are schematically illustrated in Fig. 3. The 
figure shows “snapshots” of the current sheet at several 
different times. The branch comes into existence at some 
time i c ; the points BA and BT are assumed to initially be 
separated by the distance L. The point BA moves with ve- 
locity v ba = |v ba |x. The point BT has two velocity com- 
ponents, Vbt = |vht| sin Ox + | Vbt | cos Oy, because the 
flux tubes that wrap around this point (see frames 3 and 
4 of Fig. 2) exert both axial and transverse components 
of magnetic pressure, driving the point both forward and 
downward. The branching phase is assumed to terminate 


when the point BT reaches the cathode (indicated by time 
tf in the figure). The predicted terminal canting angle 
(Of) is calculated using the axial separation of the points 
BA and BT (s in the figure) at time tf. For a given elec- 
trode separation distance h: 

Of = tan -1 , (1) 

rU 

s= (|v ba | - |v bt | sin0)dr , (2) 

Jo 

rU 

h — L = / | Vbt | cosfldr . (3) 

Jo 

Assuming |vba| and |vbt| are known, Eqns. 1-3 can be 
numerically integrated to yield the predicted canting an- 
gle. The values of |vba| and |vbt| must come from phys- 
ical models that describe current sheet propagation. In 
the next section we make arguments for the appropriate 
choice of values for these constants. 
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Figure 3: Schematic of current sheet canting angle calculation elements (rectangular geometry). 


4 Current Sheet Propagation 

4.1 Propagation Speed 

An ideal current sheet forms at the discontinuity between 
a vacuum region filled with expanding magnetic flux and 
a conductive, relatively flux-free plasma. The extent to 
which the magnetic field penetrates into the plasma de- 
pends on the thermodynamic state of the plasma. Clas- 
sically, two extreme cases can be distinguished. If the 
plasma is highly conductive, the penetration of the fields 
is confined to a thin layer of thickness on the order of the 
electron skin depth and the plasma is pushed by the mag- 
netic pressure (snowplowed) with a characteristic speed 
(the Alfvdn speed) [ 1 7] 


V A = ,. , (4) 

y/4iTnimi 

where B is the magnetic inductance, n* is the ion den- 
sity, and m.i is the atomic mass of the ions. If, on the 
other hand, the plasma is highly resistive, the magnetic 
field simply diffuses through the plasma (without impart- 
ing significant momentum) with characteristic speed[17] 


vr> = 


(?T) 

4 7rL D ’ 


(5) 


where c is the speed of light in vacuum, 77 is the plasma 
resistivity, and Ld is the characteristic diffusion length. 

4.2 Hall Effect Enhanced Field Penetration 

Fruchtman et a/[18, 19, 20] describes a mechanism for 
fast magnetic field penetration into a plasma which is in- 
dependent of the resistivity, which he calls the “convec- 
tive skin effect”. The mechanism results from the Hall 
electric field, which allows the magnetic field to penetrate 
the plasma with characteristic speed 

s y (6) 

47 :L C 

where — e is the electron charge and Lc is the characteris- 
tic length of a density gradient (transverse to the propaga- 
tion direction) in the plasma. The term t)h = B/n e ec can 
be thought of as a “Hall resistivity”. As with conventional 
plasma resistivity, the Hall resistivity allows the magnetic 
field to rapidly propagate through a plasma without im- 
parting significant momentum to the ions (the Hall re- 
sistivity is, however, non-dissipative). The Hall effect- 
enhanced magnetic field penetration will occur when a 
density gradient occurs in the (pushed) plasma with char- 
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acteristic dimension of the order of the ion skin depth. 


Si = 



(7) 


Also, a parameter equivalent to the magnetic Reynolds 
number can be defined, which is the ratio of the plasma 
pushing speed (v a) to the field penetration speed (vq): 



b 

\Z47t nnij 
c 2 (B/nec) 
4 nL 



( 8 ) 


Equation 8 shows us that when Lc is less than the ion skin 
depth, field penetration dominates over pushing. 


this mass motion can lead to anode starvation and the 
Hall-enhanced B field penetration described in the pre- 
vious section. The speed of the current sheet along the 
anode will be given by Eqn. 6. In order to calculate the 
Hall-convective speed, the scale length, Lc, for the dia- 
magnetic drift induced density gradient normal to the an- 
ode must be known. In Fruchtman’s analytical theory, this 
scale length is defined as Lc = {(d/dj/)ln[n e (j/)]} _1 . In 
order to make canting angle predictions either the func- 
tional form of the electron distribution near the anode 
must be provided from a model, or scale length estimates 
can be made using experimental data. 

5 Model Results 


4.3 Diamagnetic Drift and Anode Starvation 

The sequence of events, described in section 3.1, that lead 
to a canted current sheet are initiated by depletion of the 
conducting plasma in the vicinity of the anode. In this 
section we describe one mechanism that may drive this 
depletion process: diamagnetic drift. 

One of the most conspicuous physical traits of current 
sheets are the large gradients in both the field strengths 
and thermodynamic state variables that occur across it. It 
is therefore natural to first check the influence of these 
abrupt transitions on the phenomenon of interest. In cur- 
rent sheets the gasdynamic pressure and magnetic field 
gradients are generally quite large. Plasma pressure gra- 
dients transverse to magnetic fields gives rise to motion of 
the plasma in a direction perpendicular to both the B field 
and the pressure gradient. For example, a simplified ver- 
sion of the fluid momentum equation for the j t h plasma 
species gives the component of velocity perpendicular to 
the magnetic field[21]: 


E x B 



UE 


Vpj x B 

qjfijB 2 

u d 


(9) 


The first term in Eqn. 9 is the E x B drift (ug) and the 
second term is the diamagnetic drift (uj). The direction of 
ue is the same for both ions and electrons but the diamag- 
netic drift direction is seen to be charge dependent. With 
the anode on top as illustrated in Fig. 3, E is directed in 
the —y direction, and B is directed in —z direction so that 
ue is directed in the +x direction. Similarly, with Vp in 
the +x direction, for ions, u,j is directed in the — y direc- 
tion, or away from the anode throughout the current sheet. 

Since most of the plasma inertia is carried by the ions, 
the simplified analysis above indicates that current sheets 
with large pressure gradients in the axial direction will in- 
duce fluid flow away from the anode. We propose that 


Referring again to Fig. 3, if we assume that the branch- 
trunk interface motion is governed by plasma pushing 
(|vbt| = va), and that the branch-anode interface motion 
is governed by Hall-convective transport (|vb a | — v c ), 
then Eqns. 1-3 can be integrated to yield the canting an- 
gle for a given propellant atomic mass, ion number density 
and characteristic density gradient length. 

The model results for the parametric dependence of 
canting angle on propellant atomic mass and number den- 
sity are plotted in Fig. 4. The axis ranges were selected 
to encompass the mass and density regime of interest 
for pulsed plasma thrusters. We have fixed h=5 cm and 
L c =0.1 cm, which are typical values for a pulsed EM 
accelerator, as discussed below. The model predicts in- 
creased canting with increasing atomic mass and decreas- 
ing number density. 

5.1 Comparison with Experiment 

In an earlier study[5] we measured the canting angle of 
current sheets using a variety of different propellants (hy- 
drogen, deuterium, methane, helium, neon, argon, kryp- 
ton, and xenon) and pressure levels (75, 200, and 400 
mTorr). The accelerator was rectangular in geometry, 
with dimensions 5x10x60 cm (height x width x length). 
We also measured the current sheet speed, B-field, elec- 
tron density, and electron temperature. This data can be 
used to estimate va and v c . The diamagnetic drift speed 
was estimated to be |ud| ~ C?(10 -3 ) m/s. The charac- 
teristic time for the branching phase was observed to be 
U ~ C>(10 -6 ) s, which implies that diamagnetic drift is 
capable of displacing plasma away from the anode about 
|u<i|f6 ~ 0.1 cm during the branching phase. Therefore, 
we take L c ~ 0.1 cm as the order of magnitude of the 
near-anode density gradient. The magnitude of this char- 
acteristic length is perhaps not surprising, since the ex- 
perimental measurements yield a comparable estimate for 
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Figure 4: Canting angle dependence on atomic mass and ion density (h=5 cm, L c = 0.1 cm) 


the axial characteristic length, n/Vn. Using the electron 
density data we found that n/Vn ~ 0.1 cm for all cases, 
which is the same order of magnitude as our estimated 
normal characteristic length. 

Using the diamagnetic drift-induced characteristic den- 
sity gradient ( L c ~ 0.1 cm) and the average experimen- 
tally measured number density (~ 1 x 10 16 cm -3 ) we can 
calculate the predicted canting angle as a function of the 
propellant atomic mass. This curve, along with all of the 
experimental datum points are shown in Fig. 5. The fig- 
ure shows that the model captures the general trend of the 
experimental data; the model predicts that the canting an- 
gle will initially rapidly increase with propellant atomic 
mass and then taper off toward a more asymptotic trend 
for higher atomic mass propellants. These are the same 
trends that are seen in the experimental data. 

It is not surprising that the model does not give pre- 
cise quantitative agreement with the data for the higher 
atomic mass propellants. The model predicts the “worst 


case scenario”, or maximum expected canting angle, be- 
cause the simplification in the model that allows the point 
BA to propagate with velocity v c along its entire trajec- 
tory probably does not completely capture the complex- 
ity of the actual near-anode behavior. In the real accel- 
erator, the branch anode attachment is likely to transition 
continuously from the Hall-convective behavior back to 
plasma pushing - as it picks up speed and the incoming 
propellant obliterates the anode density gradient. If we 
were to modify the model to divide the branch trajectory 
into, say, two parts (half Hall-convective and half plasma 
pushing), the effect would be to cause the model curve to 
bend-over more rapidly and come into closer agreement 
with the data. But, in light of the rather crude order of 
magnitude estimates used throughout this analysis, fur- 
ther refinement of the model is not warranted. The im- 
portant point is that the simplest model, which embodies 
the essence of the proposed physical processes, accurately 
picks up the trends in the experimental data and gives 
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propellant atomic mass [amu] 


Figure 5: Comparison of model and measured current sheet canting angle[5] versus propellant atomic mass. 


fairly close quantitative agreement. 

6 Discussion and Conclusions 

The proposed current sheet canting model can be summa- 
rized as follows. Physically, Fruchtman’s model shows 
that the Hall effect can lead to an effective mechanism 
for magnetic field transport when a non-uniformity (i.e., a 
density gradient) appears in a plasma that is being pushed 
by the field. This argument is formally derived through 
the fluid equations. In particular, inclusion of the Hall 
term in the generalized Ohm’s law leads to a term in the 
equation for the time-evolution of the magnetic field that 
allows for the rapid convection of magnetic field through 
the plasma, without displacement of the ions, with char- 
acteristic speed v c . This competes with the convection 
of the magnetic field where the ions are displaced, that is 
plasma pushing, which has characteristic speed v\. Con- 
sideration of the relative magnitude of each effect reveals 
that the Hall-convective penetration becomes comparable 
in magnitude to plasma pushing when a density gradient 
of scale length comparable to the ion skin depth appears. 

In our fast field penetration canting model we propose 
that such a density gradient does arise along the anode of 


the accelerator due to diamagnetic drift, allowing local- 
ized rapid penetration of the magnetic field and the for- 
mation of two distinct current carrying structures, which 
we call the branch and the trunk. Current sheet canting 
results from a disparity in the axial propagation speed be- 
tween the branch-anode and branch-trunk interfaces. The 
branch-anode interface moves faster and the current chan- 
nel connecting the anode to the trunk forms a canted cur- 
rent sheet. The canting angle continuously evolves from 
zero to its terminal value during the branching phase. The 
severity of the canting depends on the relative magnitude 
the branch-anode and branch-trunk interface velocities. 
The variations in the ion skin depths among the d ifferen t 
propellants (the ion skin depth is proportional to ^/ini/n,) 
leads to different Hall-convective field penetration speeds, 
and hence different final canting angles, as illustrated in 
Fig. 5. 

This study’s practical recommendation is that future 
GFPPTs should use low atomic mass propellants at high 
pressure to avoid performance losses due to current sheet 
canting. Hydrogen would seem to be the best choice; 
however, hydrogen is not an ideal PPT propellant, due 
to the difficulty of handling cryogenic propellants on a 
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spacecraft. To practically exploit the benefits of low cur- 
rent sheet canting which results from the use of hydro- 
gen at high pressure, alkanes, with their hydrogen-rich 
CnH-in+i structure, may be a natural choice. Methane 
was tested as part of the experimental study and was found 
to have the same reduced canting behavior (at higher pres- 
sures) as hydrogen. Further tests are needed on longer- 
chain hydrocarbons. If butane, for example, is found to 
exhibit similar behavior, it will be an appealing GFPPT 
propellant on two levels. In addition to the aforemen- 
tioned canting benefits, butane can be stored as a liquid 
under relatively low pressure, at room temperature; there- 
fore, a butane propellant system will have a much smaller 
specific volume (smaller fuel tank and feed system) than 
a high pressure gas system, and provide a lower overall 
system a. 
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